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Stochastic mean-field equations are derived applying the time- dependent 
variational principle to a quantum many-body system in interaction with 
a classical heat bath. This approach is tested on the case of a charged par- 
ticle in thermal radiation field, and is used to include frictional and random 
forces in the TDHFB equations. Dynamical estimates of the diffusion coef- 
ficient for the quadrupole deformation are obtained. 
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1 Introduction 



The simulation of the non-equihbrium fluctuation phenomena observed in 

hot nuclei represents a fundamental problem of the many-body theory, de- 
manding stochastic extensions of the kinetic or mean-field equations. 

At the level of the semiclassical description provided by the phase-space 
density, the fluctuations can be accounted by the Boltzmann-Langevin ap- 
proach [1]. Within this approximation the nuclear Boltzmann equation is 
modified, so that a random component (Brownian force [2]) is included in 
the collision term. 

The time-dependent variational calculations on coherent states trial man- 
ifolds, as time-dependent Hartree-Fock (TDHF), or antisymmetrized molecu- 
lar - dynamics (AMD), lead to deterministic mean- field equations, but recipes 
to introduce a random behavior have been proposed [3, 4, 5, 6, 7]. The ba- 
sic assumption is that some decoherence mechanism exists, interrupting the 
deterministic evolution of the quantum wave packet by spontaneous collapse 
to the neighbouring states of the same trial manifold. The collapse rate is 
obtained using the Fermi's golden rule, for a suitable residual interaction. 

The mechanism of the random spontaneous collapse induced by environ- 
mental (residual) fiuctuations remains obscure, beyond the framework of the 
quantum many-body theory. This ambiguous mechanism can be avoided 
if the residual interactions and the unretained degrees of freedom are ac- 
counted from the very beginning in the time-dependent variational principle. 
Particularly suited for this purpose is the model of bilinear coupling between 
a particle and a heat bath of harmonic oscillators, applied with success to 
derive the classical Langevin equation [8] . This model can be used also to in- 
clude stochastic forces in the mean-field equations, and was apphed to study 
the effects of pairing, dissipation and temperature on the giant quadrupole 
resonance (GQR) [9]. 

The dynamical approach to the coupling between a quantum many-body 
system and the heat bath is presented in Sec. 2, together with a test on 
the important case provided by a charged particle in blackbody radiation 
field. The stochastic mean-field equations are used in Sec. 3 to calculate the 
diffusion coefficient of the quadrupole deformation coordinate. Conclusions 
are drawn in Sec. 4. 
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2 Brownian quantum dynamics 



The phenomenology of the Brownian diffusion in classical mechanics can be 

described by two external force terms in the Hamilton equations of motion 
for the particle: a random force with zero mean, ^(t), (the noise), and a 
dissipative friction force, /. The type of friction is determined by the memory 
function T{t), or the spectral density of the environment, J{uj), so that for a 
linear dissipation mechanism 

m = - /' r(t - . m - - r d.£M£5^fM , (1) 

Jo TT Jo UJ 

with Q the particle velocity. The memory function gives also the noise cor- 
relation function, and ^{t) should obey the fluctuation-dissipation theorem, 

«mm»-kBTr{t-t') , (2) 

where << .. >> denotes the average over an ensemble of trajectories. 

Noise and friction forces having these properties may be derived follow- 
ing a Hamiltonian treatment of the whole system consisting of the particle 
in interaction with a heat bath of Nc harmonic oscillators via the bilinear 
coupling term [8] 

Hcoup — Q 5Z ^j^j 1 (^) 

i=i 

where Cj are coupling constants and qj the time-dependent bath coordinates. 

The advantage of the bilinear coupling model over the phenomenological 
Langevin approach becomes striking when the particle is a quantum ob- 
ject. For example, below a certain crossover temperature, the mechanism of 
chemical reactions changes from thermal activation to tunneling, becoming 
sensitive to destruction of the phase coherence and dissipation produced by 
environment. The temperature effect on the tunneling rate can be accounted 
by the thermal average of the quantum penetrability factors [10]. However, 
dissipation requires a more elaborate treatment, based essentially on Eu- 
clidean path integral calculations for the bilinear coupling Hamiltonian [11]. 

The problem of dynamical coupling between a quantum system and its 
environment is even more complex, having not yet a satisfactory solution. 
Thus, no standard method exists to treat a mixed system composed of a 
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quantum particle interacting with a classical heat bath. In quantum gravity 
the similar problem of a quantum system coupled to the classical space- 
time foam could be treated by Hamilton-Heisenberg equations derived from 
a variational principle [12]. In the quantum many-body theory, the collective 
mean-field obtained from time-dependent variational calculations has already 
a hybrid character, including both quantum and classical aspects. Thus, it is 
natural to define its dynamical coupling to a classical heat bath of oscillators 
using the variational equation 

5 / dt[Y,miqi)' + Wm-i'ao + iibm]^o (4) 

1=1 

where Hq is the Hamiltonian operator for the isolated quantum system, and 
H5 contains the bath energy plus the coupling interaction. 

Let us consider variations of the quantum wave function ^ in Eq. (4) 
restricted to a finite dimensional trial manifold S = parameterized 
by 2N coordinates X = {x^}, i = 1,2N, so that the matrix uj^ = [cu,^(*)], 
defined by 

u;f.{^^^) = 2hlm{di'i/\dj'i/) (5) 

is non-singular [13]. Supposing a coupling interaction Hcoup — '^J2f=i CjQj, 
linear in the coordinate operator K, the equations of motion determined by 
variations in Eq. (4) with respect to X{t) and the bath trajectories qi{t) are 

for the quantum system and 

Qi^— , Pi^-miu;^qi-Q{^\K\^) (7) 
rrii 

for the classical oscillators. The classical equations can be solved in terms 
of the unknown function of time = (\E'|K|\E'), and when their retarded 
solution is inserted in Eq. (6), a Langevin-like equation is obtained. 
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Here is given by Eq. (1) with instead of Q. The term 
K[(5,/,(0)r(t) — (5^(t)r(0)] wiU be neglected, supposing that it can be ac- 
counted by renormahzing Hq. 

If S is the whole Hilbert space of the quantum system, then Eq. (8) 
becomes a non-linear time-dependent Schrodinger equation 

^h^ = [u,-\^m + U)]'^ ■ (9) 

The non-linearity appears in the friction term Z^, accounting for the back- 
reaction of the environment on the quantum system. When friction is ne- 
glected, the rate of the noise-induced excitation between two eigenstates \Ei) 
and \Ef) of H can be easily calculated in the first-order approximation for 
the time-dependent perturbation represented by ^noise = The re- 

sult is 



|2 

■liiTif 



h f i{s)e'^n'ds\^ (10) 
t Jo 



and can be further averaged over the statistical ensemble before the time- 
integral making use of Eq. (2), so that 



2 |2, j.J{^fi) 
-2 -K/i /tfii 

^ UJfi 



rf''^=«rfi»^^\Kf,fkBT^-^ . (11) 



Here K/j = {Ef\K\Ei) is the matrix element of the couphng operator K, and 
Ufi = \Ef - Ei\/h. 

The rate provided by Eq. (11) may be used to test the formalism pre- 
sented above on the case of a single charged particle in thermal radiation field, 
considered as a classical heat bath of oscillators continuously distributed in 
frequency. In the dipolar approximation, the coupling of the electromagnetic 
field tojhe charge e has the form = -d E = -d^E^ - dyEy - d^E^ 

where d = er is the electric dipole operator, and E the radiation electric field. 
For isotropic radiation H^f^g is a sum of three terms of the form — K,^ as- 
sumed above, each with the spectral density J™'^{u) = 2he^u}^{n)i_o/ {3c^kBT). 
The total transition rate Ty^, is the sum of three terms provided by Eq. (11), 
and the result 

r? = ^|?/.lH(^W' {n). = l/{e'-/'-^-l), (12) 



5 



coincides with the known transition rate in thermal radiation field [14]. At 
high temperatures the memory function corresponding to J^"''^{u) can be 
expressed in terms of the second derivative of the delta function, r{t) — 
— 4e^5(i)/(3c^), and the related friction force 



frad _ /-iqN 



is the classical radiation reaction in vacuum [15]. 

In the case of Ohmic dissipation r{t) — 2^^^^), with the static 
friction coefficient, and Eq. (8) becomes 



dx'^ dt 

3 TDHFB-Langevin dynamics and the shape 
diffusion coefficient 

An important application of the stochastic mean-field techniques concerns the 
calculation of the shape diffusion coefficient. The thermal equilibrium state of 
a highly excited nucleus can be represented [7] by an ensemble of incoherent 
self-consistent Hartre-Fock configurations |^a;), with energies Ek — (Ho)*^. 
and quadrupole deformations 



Within this ensemble the transitions induced by the residual interactions 
occur with a rate 

2n 

T/i = -^\{^f\Hres\^,)\'S{Ef - E,) (16) 

and a diffusion coefficient of the quadrupole deformation can be defined by 
[7] 

= - A)'r^. , (17) 
/ 
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where the sum involves states |^/) different by in the occupation number 
of two orbital levels. However, the value predicted by this calculation is only 
a fraction {^i 10%) of the estimate provided by the Einstein formula^ 

= ^ (18) 

with the friction coefficient 7^ extracted from experiment. 

A dynamical calculation of Dj3 can be performed using an ensemble of 
Brownian trajectories P{t) generated by Eq. (14) for a nucleon system in 
contact with a heat bath. The study of the axially symmetric quadrupole 
shape vibrations in the sd shell nucleus ^^Si including pairing, friction and 
temperature shows that the occupation numbers p,.,r = 1,6 of the orbital 
levels depend on time and may switch between different configurations, but 
each configuration is quasi-stable for relatively long periods, despite the shape 
fluctuations. For ^^Si this switching means a change of {pr} from the oblate 
configuration of the metastable ground state (mgs), to an unstable 

prolate configuration, {p*}. The rate of these oblate^prolate configuration 
transitions can be defined in terms of the mean-first-passage-time [8], r*, 
which is the average time elapsed until the system switches the configura- 
tion for the first time. Therefore, a diffusion coefficient of the (quasi-stable) 
quadrupole deformation can be defined by the average 



1 (fjb _ iDb\2 

- ]^ E J' > P^^P\0) , P'f^P\r;) (19) 



where the sum is taken over trajectories generated with random initial 
conditions at a fixed energy Ei. For such calculations, the trial manifold 
S in Eq. (4), which is the best suited to account for the occupation num- 
ber degrees of freedom, consists of a set S^^^ of normahzed Hartree-Fock- 
Bogolyubov functions. In this case, without environment coupling, Eq. (14) 
corresponds to the deterministic time-dependent Hartree-Fock-Bogolyubov 
(TDHFB) equations. The sd shell contains Ng = 12 states, and S^^^ is 
parameterized by Ns{Ns — 1) = 132 variables {a;*}. This number is large, 
and therefore it is convenient to restrict S^^^ to a submanifold with lower 
dimensionality, accounting only for the quadrupole and occupation numbers 



^When xt is a classical Brownian trajectory, limt^oo « (xt—xo)^ » /t = 2kBT/jx 
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degrees of freedom. For the 6 protons and 6 neutrons of ^^Si outside the ^^O 
core, the 13 2- dimensional manifold S^^^ can be reduced by the cranking 
technique [17] to a 14-dimensional submanifold corresponding to 7 degrees 
of freedom: (3 and the occupation numbers = 1; 6 of the sd shell orbital 
levels. 

The model Hamiltonian Hq includes the spherical harmonic oscillator 
one-body term hg and the quadrupole-quadrupole plus pairing two-body in- 
teractions: 

Ho = ho-xK'-^P^P ■ (20) 

Here K = Qq/cq, Cq = {b/ATrY^^h/muJo = 26.17MeV-fm7/ia;o, coq = AlA-^/^ 
MeV/h, and the interaction strengths are chosen to be % = 0.186 MeV and 
Gpair — 1-23 MeV. The mgs configuration corresponds to an oblate 

shape = —0.43), and is obtained by frictional cooling [18]. This procedure 
ensures the self-consistency of the quadrupole and pairing mean-fields [9]. 

The operator K is also the coordinate appearing in the coupling to the 
heat bath. The corresponding friction coefficient is related to by jk = 
0.00837/3, and at T = the rate F of the energy dissipation is the same as 
the decay width of the GQR, Fgqr ~ 4.5 MeV, when ~ 6Ah. 

The average in Eq. (19) was calculated using 300 trajectories starting at 
the GQR energy of 19.1 MeV, with random initial conditions. When ksT is 
between 1.5 and 2.5 MeV, and jp is near GAh, then can be accurately 
interpolated by the formula 

D*(T) = ^^'^ (21) 

^ u + vkBT + -fp{a + bkBT) ^ ' 

where u = 0.27h, v = l.lh/MeV, a = 0.14 and b = O.lMeV'^ This result 
indicates that the effective friction coefficient for shape diffusion with config- 
uration change, 7^ = m -|- vksT + ^^{a + hkBT) is not the same as the fixed 
model parameter 7^3, but depends on temperature. For ksT — 2.5 MeV, 
Eq. (21) predicts a diffusion coefficient — 453 keV/h, about a factor 
ten larger than the microscopic estimate [7] of 30 keV/h in the neighboring 
nucleus ^^Mg. Though, assuming an dependence of on the mass 
number, for ^^^Er one obtains ~ 80 keV/h, which is close to the value of 
~ 70 keV/h extracted from experimental data. 
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4 Conclusions 



The variational principle applied to a quantum many-body system and its 
thermal environment provides stochastic mean-field equations (Eq. (8)). De- 
pending on the trial manifold, a non-linear Schrodinger-Langevin equation 
(Eq. (9)), or a system of TDHFB-Langevin equations can be obtained. The 
example of a charged particle in blackbody radiation field shows that this 
formalism predicts correctly the transition rate and radiation reaction. 

The stochastic mean-field equations (TDHFB-Langevin) have been ap- 
plied to calculate the diffusion coefficient of the quadrupole deformation (Eq. 
(21)). The effective friction coefficient 7^ for shape diffusion with change of 
the internal configuration is not the same as the model parameter 7^ and de- 
pends on temperature. The results of the present calculations within the sd 
shell are consistent with the measured values of the GQR decay width and 
shape diffusion coefficient, proving the relevance of the TDHFB-Langevin 
formalism. 
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